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1. GENERAL OUTLINE OF THE MODEL 


Although the core of the model used in 
"Scope for Growth" is of the classical input-output type, 
some adaptations have been introduced to enable it to be 
used as a tool for prospecting the future. We will begin 
by describing the model in broad outline, using the 
familiar terminology of input-output models. Next we will 
focus on the adaptations. In Section 2 the complete model 
is presented in detail. 


In traditional input-output models, the 
following equations can be defined for each sector i of 
the economy: 


X; = = X; jt C; +I; + E; 

Where gross output of sector i (=X; ) equals: 

- the sum of intermediate outputs of sector i to be used 
in the same sector or in other sectors (sum over j of 
Xij)? É 

- plus the output of sector i that is consumed as end 
products (C;), 

- plus the output of sector i to be exported (Ej): 

- plus the output of sector i to be used as capital 
goods in other sectors (I;) 


These equations are known as the Leontief equations. 


Introducing the technical coefficients a; 33 


ae S qo 7s Mae a 
representing the amount of input from sector i needed to 
produce one unit of output in sector j and substituting it 
in the equations above, we get (in matrix notation): 


(1) X = A X + C + E + I 
Xi 
where X is a vector of the form X 
Xn 


The vectors C,E and I have the same form. A is a matrix 
with characteristic elements a; j" 

Usually, assumptions are made for changes in 
one of the demand components and the corresponding (new) 
levels of sectoral production needed to satisfy the 
demand, are then calculated. This is done by rewriting 
equation (1) as follows: 


X = (ID-A) Í . (C + E + I) 


where ID = identity matrix. 


The model can be made dynamic by linking the 
levels of investments in one period of time to the 
production capacity in the following periods. The needs 
for imports are usually defined in an extra set of 
equations, and are related to levels of production, 
consumption and investment. Also, labour equations are 
usually introduced, linking production levels to 
employment. 


Special features 


Our model, TEB-4 (Technical-Economic Bookkeeping, 
version 4) differs in three main respects from classical 


input-output analysis. 


l. Instead of varying any of the demand components in 
order to see what happens to the production level and 
its sectoral composition, in TEB-4 upper and lower 
limits are chosen and set for consumption, exports and 
indirectly for investment of goods and services in the 
different sectors. These limits reflect expectations, 
wishes or technical necessities, as explained in 


“Scope for Growth". Consequently, there will not be 
only one result for the level and composition of 
production given by the model, but a large number of 
production levels and production compositions will 
satisfy the demands. Thus room is created for 
optimization. 


2. In TEB-4 the Leontief and other equations have an 
additional index t for each of the prospected years in 
the future. As we prospect the next ten years, this 
increases the number of equations in the model by a 
factor of 10. | 


3. The third new feature of our model concerns the 
environmental equations, These equations are not a 
traditional part of input-output analysis but are in 
keeping with it, and relate the emission of pollutants 
on the one hand to production and consumption levels 
on the other. Besides the 18 traditional sectors that 
are defined in the model (see Annex 1), 5 more sectors 
are defined where pollution abatement activities are 
undertaken. These sectors do not produce output like 
the other sectors (the X variables in the model), but 
produce abated pollution value. 


In matrix notation, the following equations form the main 
framework of the model used: 


(1) Leontief 


X = A. X + C + E + IÏ 
(2) Investment 
I = D. X + K.AW 


(3) Imports 
M = M .X + M . C+M . I 


(4) Labour 


(5) Pollution 

where X, C, E, I, AW, M, L, V are vectors of model- 
variables, A, D, K, Mz’ Mo’ M. Vy Me are 

matrices with the technical, depreciation, capital, import 
and emission coefficients, and gE anā AP are 
parameters for the embodied and the disembodied technical 
progress. j 

These equations will be presented in detail below. 


2 FORMAL PRESENTATION OF THE MODEL 
2.1 Introduction 


A dynamic input-output model has been 
formulated, capable of describing the changes in eighteen 
traditional economic sectors and five (additional) 
pollution abatement sectors over a future period of ten 
years. In essence the input-output framework ensures that 
two kinds of requirements are met. First there is the 
requirement that everything produced in a sector is also 
sold. The sales possibilities include supplies to other 
sectors or within the sector itself, deliveries of capital 
goods, deliveries to consumers, or exports. No account has 
been taken of changes in stocks. Secondly there is the 
requirement that, according to the technical relationships 
(the production functions) laid down in advance, the 
factors required to produce that output are available. 
These consist of imports of raw materials and 
semi-manufactures, inter-industry sales and labour and 
capital. The inter-industry sales’ are both input and 
output. In a dynamic input-output model of this kind, the 
same applies to domestically produced investment goods. 
Capital goods supplied in a given year are assumed to 
provide capital services for the recipient sector in the 
following year. Thus investment per sector has been 
related to the required level of productive capacity. 

If the final demand (i.e. consumption, 
exports and deliveries of capital goods) per sector is 
given, the required level of labour, imports, investment 
and inter-industry sales can then be deduced. The 
calculation of these requirements is a traditional 
application of input-output analysis. No fixed values for 
consumption and exports have, however, been selected in 
this model; instead intervals have been specified within 
which they can move in time, thus providing room to 
examine potential future developments by means of 
optimizations. With respect to the domestic and foreign 
demand for end products, margins have been introduced by 
setting upper and lower limits to the sectoral growth of 
consumption and exports. In this way, uncertainties with 


respect to the behaviour of consumers and exporters are 
reflected in the spread of these margins. 

The model is supplemented by equations 
describing the emission of selected environmental 
pollutants as a function of sectoral production and 
consumption levels. The requirements in terms of inputs 
(labour, capital, imports, raw materials) to control 
pollution are specified. Whether or not these additional 
abatement activities will be executed, however, depends on 
the weights given to the various goals in the optimization 
process. Finally, a number of goal variables have been 
included, ranging from economic, socio-economic to 
environmental objectives. 

“Though the structure of the model is 
relatively simple, its size is considerable because 
23 economic sectors have been distinguished and we want to 
prospect the next 10 years. It comprises 2213 equations 
-("rows" in linear programming terminology), approx. 1400 
variables (“columns”) and approx. 12,500 non-zero 
coefficients. All equations are linear, all variables 
non-negative. 792 of the equations are in the form "equal 
to", the others being inequalities. In the following 
sections a formal presentation of the model is given. 
Eight blocks will be distinguished: (i) gross output 
(ii) investment and production capacity, (iii) exports, 
(iv) imports and the balance of trade, (v) consumption, 
(vi) labour, (vii) the emission and possible abatement of 
pollutants and (viii) the goal variables. Each block 
includes definitions and, in most cases, limits on 
variables or their changeability. 


With respect to the notation, the following 

conventions are observed: . 

- Capitals: goal variables 

- Lower case Latin letters (except a, k, d and q): other 
model variables 

- Lower case Greek letters and a, k, d and q: 
coefficients 

- Index i: sector of origin (row in the input-output 
table); the terms "(economic) sector" and "industry" 
are used interchangeably 


(1). 
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- Index j: sector of destination (column in the 
input-output table) 
- Index t: year; t=o represents the base year 
- Index k: type of pollutant 
- Upper index*: exogeneously determined value 


- A z annual change 


All variables are in million guilders and 


constant prices, unless otherwise stated. All coefficients 


are in unitary dimensions, unless otherwise stated. Under 


the formulas presented in this section a short description 


of the variables and coefficients is given. For the 
extended definitions see Annezes 4 and 5. 
2.2 Gross output 

According to the Leontief definition 


equations, the gross output ("production") of any industry 


must add up to the sales of intermediate deliveries'to 


by households and government, plus exports, plus sales 
investment goods. This holds true for each of the 18 


. itself or to other industries, plus domestic consumption 


“of 


conventional economic sectors for any year, bringing the 


total number of equations to 180 (for the sectoral 
breakdown see Annex 1). 


He 
i 


With as variables: 


1,...18 
1,...10 


x = gross output, domestically produced 

c = consumption by households and government, 
domestically produced 

e, = exports 

r = interest margins of banks (not equal to zero 


only for the Business Services sector). 
gross investment, domestically produced 


He 
ti 


and coefficients: 

ai,j = technical coefficient: 
needed to produce one unit of output of 
industry j ` 


sales of industry i 


(la) 


243 nv ment an d i a i 

Gross domestically produced investment is 
divided into replacement investment (its allocation being 
determined by the depreciation coefficient matrix) and 
expanding investment (its allocation being determined by the 
capital coefficient matrix). Furthermore, for some recipient 
industries investment programmes are predetermined, 
resulting in exogenous variables for the investment goods 
producing industries and, as will be shown later, in 
exogenous variables in the import equations. Predetermined 
investment programmes are formulated for the sectors Mining 
& gas distribution, Electricity & water, Housing, and Public 
administration & defence. For these recipient industries the 
description of the allocation mechanism by capital and 
depreciation coefficents is dispensed with: in other words, 


for j = 2,10,15,18, Kise = dijt = 0. 
= S je i= 1 18 
s en ` ` j l= eee 
ipes B Bae Peke Me eRe a 
j=1 ` l j=1 t= 1,...10 

with as variables: 

Wa = production capacity 

io = exogenous investments 

and coefficients: 

dij = depreciation coefficient: sales of investment 
goods produced by industry i needed to 
maintain one unit of production capacity in 
industry j 

ki,j = capital coefficient: sales of investment goods 
produced in industry i needed to expand the 
production capacity in industry j by one unit. 
Although time indices are attached to the 

technical, depreciation and capital coefficients - implying 

that changes in time can be allowed for if such information 

is available - in the actual computations the coefficients 

are generally assumed to remain constant over the next ten 

years. 


Investments appear in the Leontief equations 
as a demand category; they are also relevant in relation to 
the production capacity as equations (la) show. By rewriting 


(2) 


(3) 


equations (la) with only w, on the left-hand side it is 


easy to see that capacity in year t is a function of the 


‘capacity in the preceding year, the expansion investment and 


the replacement investment. In the computer model 
equations (la) are substituted in equations (1). For the 
domestically produced part (for the required imports see 
section 2.5) the production capacity of each industry is in 
this way implicitly defined. 

We now state that production in any year may 
not exceed the capacity at the end of the previous year. 


i= 1,...23 
t= 1,...10 


To avoid too high a degree of idle capácity a lower limit of 
production in relation to capacity is set. This also 
provides against "ornamental investment" which can appear as 
an outcome if certain conflicting requirements are imposed 


‘(in earlier exercises, for instance, a high level of 


employment and low use of energy (1)). 


Ï X. “> & Ïw ” i=l,...23 


with the coefficients: ; 
= minimum average degree of utilization of capital 
(= o for i = 2, 10, 15, 18) 


Expansion of capacity cannot be unlimited, one reason 
being that there is not an unlimited supply of qualified 
labour to man the new capacity. Therefore 


° 


(4) w; | $S lta.) w i= 1,...18 i#2,10,15,18 
cT, eee ld 


with as coefficients: - 
w = maximum permitted annual expansion of capacity. 


For those industries where exogenous investment programmes 
have been specified, the future capacity is also fixed. 


(5) w. = wF i= 2,10,15,18 
t= lyo LO. 


Sectors 2, 10, 15, 18 are Mining & gas distribution, 
Electricity & water, Housing, and Public administration & 
defence. In the case of the pollution abatement sectors 
(nos. 19-23) an absolute (that is: not related to the 
capacity in the previous year) limit is set. This 
represents a financial ceiling to additional environmental 
investment expenditure.. j 


19: 7 < 6623 
= 1,...10 


i 


(6) Aw, _ < w. i 


ct 
1 


The contraction of capacity is also subjected to a limit. 
No deliberate dismantlement of capacity is allowed, 
implying for each industry a rate of fall not exceeding 
the rate of depreciation. In other words, negative gross 
investment is excluded. 


18 
MK : . , 
{( l + Tie) E Kjt] aw, tÉ j=l, ...23 j#2,10,15,18 
t=1,...10 
18 
MD 
(( 1 + aje) z S Xit 
i=l 
where: 


= import coefficient for expansion investment 
import coefficient for replacement investment. 


<2 <2 
It 


For those industries with pre-determined investment 
programmes, restriction (7) is superfluous. For the 

: „j = 0, 

restriction (7) implies a non-decreasing capacity during 
the projection period. 


pollution abatement sectors, where d; 


A final restriction concerns investment in 
the last year of the prospected period (t=10). In that 
year there is no future perspective serving as a guide for 
adjustment in production capacity. In order to avoid a 
sudden decrease in production capacity in the last year of 
the prospected period, an extra restriction is added: the 
sum of expansion investment, domestically produced and 
imported, in all sectors should at least equal the sum of 
these investments in the preceding year. 


18 18 

E {1 +dg,) Ek. ., tow, 
j=1 , t iel L,j,t 3,10 

18 18 

MK 

Ï (( 1 + q. .) k. . _ }Aw 

j=1 j,t iel i,j,t j,9 
2.4 - Exports 


Exports is one of the demand categories in 
the Leontief equations (1). As stated above, no 
traditional explanatory export equations are included in 
the model. Instead upper and lower limits for the annual 
change in the volume of exports per sector are formulated. 
They are expressed as a percentage of the exports in the 
base year of the survey (indicated by subscript t = 0). 


- 12 - 


1,...18 ix 12 


G 
bt 
li 


ĉe , t == i,t i,O 
t= 1,...10 
L isl,...18 if 12 
A e aie: 0 | j 
t= 1,...10 
with coefficients: 
u 
Ê = upper limit for export change 
ee lower limit for export change 


Some industries do not produce export products (Housing, 
Public administration & defence, and the pollution 
abatement sectors). AS a consequence the upper limit and 
lower limit are set to zero or the variable is 
non-existent (in the case of the pollution. abatement 
sector). Other sectors, like Electricity & water, produce 
in principle exportable products, but are assumed to be | 
non-exporting in the projection period. Still other 
sectors do export but have a more or less fixed export 
pattern, e.g. the exports of natural gas which are based 
on long-term contracts. In this case this information is 
incorporated in the model by specifying the predicted 


te 
course (resulting in $= 5. ) 


The export volume of the Trade sector (meaning here not 
the value of the goods but the trading margins) is equal 
to a fixed ratio of the total exports of goods by all the 
goods producing industries. 


With the coefficient: 
q, = trade margin for export of goods 


As intermediate variable the exports of all sectors each 
year is defined. 


18 
(12) (aux.) e, = x è j t= 1,...10 
i=] 
2.5 Imports and the balance of trade 


Imports per sector and for each year are 
related to sectoral production, consumption and 
investment, since each of these have an imported 
component. Investments are split up into expansion 
investment and replacement investment, with their own 
import coefficients. The index j is used to denote the 
sector for which the imports are destined (just as the 
index i is used in the export equations, relating to the 
sectors where exports originated). 


18 
MI MD i 4 
13 m. =q. Xi + : j= ee 
( ) jet Tst Jet (de = ds ze) St ë 3 Ài a 
` = 1,...10 
18 
MK MC 
R L 
Se, ae AW a At Jet i 
m, 
j,t 


with the variables: 

m = ‘imports of goods and services 

m* = exogenous imports for predetermined investment in 
sectors 2, 10, 15 and 18. 


(14) 


(15) 


- 14 - 


and, the coefficients: 

= import coefficient for 
semi-manufactures 
import coefficient for 


= import coefficient for 


import coefficient for 
non-competing consumer 


9 
9, 


me 


4 


< 


MD 
MK 


raw materials and 


replacement investment 
expansion investment 


final consumption (excl. 
goods and services. These 


relate to products which cannot be produced in the 


Netherlands). 


In order to arrive at the total volume of imports each year, 


the imports of non-competing consumer goods and services 


must be added. A fixed ratio between these goods and 


services and the total domestically produced consumption is 


assumed. 
23 18 
(aux.) m. = Ï m. gen a iN Ï c è 
jar d'E, jie) 


with the variable: 
M= 


and the coefficient: 
MN 


Q 
A . 
and services 


total imports in year t 


t= 1,...10 


= import coefficient for non-competing consumer goods 


The auxiliary equations (12) and (14) combined give us the 


annual balance of trade. Annual upper and lower limits are 


placed on this balance, as well 


surplus over time. 


U 

= < $ 
e, 7 Mes "ce 
L 

-m > 6 
e, t * `t 


as to the average deficit or 


(17) 


(18) 


(19) 


(20) 


1 U 
— - < 
10 ae m.) S. ó 
‚to . 
— I (e - m ) > ô” 
10.24, t z 


with the coefficients: | 
bre upper limit for annual balance of trade 
lower limit for annual balance of trade 


w 
li 


en 
É 
u 


upper limit for average value of the balance of 


trade 
S ë lower limit for average value of the balance of 
trade 
2.6 Consumption 


Like the sectoral exports, the domestically 
produced consumption per sector is one of the demand 
components of gross output as defined in equations (1). And 
here too no behavioural equations are introduced into the 
model. Instead margins for the annual change in sectoral 
consumption are specified. These changes are expressed as a 
percentage of consumption in the base year (c; o)’ 


Ac, < yY 


it = 'i,t “i,O = 1,...18 i#l2 


Tsss LO 


d H 
i 


L a j 
bC: 2 to .i= 1,...18 API 
t= 1,...10 
with the coefficients: 
a 
Y = upper limit for consumption change 
É = lower limit for consumption change 


Like exports, the consumption of trade services is assumed 
to form a fixed ratio of the total domestically produced and 
imported consumption of goods: 


(21) 


=- LU - 


with the coefficient: 
oT 


a, = trade margin for consumer goods 
2.79 Labour 


The demand for labour is a function of 
production in the individual sectors. On account of 
differences in the level and trend of labour productivity, 
the relation varies per sector. It may, however, be assumed 
that the productivity of labour does not decline in any 
sector. In determining the demand for labour, account is 
taken of the fact that more is invested if production ina 
given sector increases rapidly thereby resulting in a 
rejuvenation of the capital goods and hence in higher labour 
productivity. To model this phenomenon it is necessary to 
distinguish between two kinds of labour-saving technical 
progress. The impact of disembodied technical progress, 
resulting from more efficient organisation, a better 
educated labour force, learning effects etc., is independent 
of the rate of growth in production. The labour-saving 
effect of technical progress embodied in equipment on the 
other hand, depends on the average age of the capital goods 
installed and is therefore influenced by the previous rate- 
of expansion of the particular industry. It can be shown 
(see Annex 2) that by assuming a survival function for 
equipment and by using a vintage approach the following 
approximate formula can be deduced. 


(22) 


(23) 


L = : 
i,t D i,tel 
LO Ae t= 1,...10 
ee... BGR (L 1114 3É 249) 
l - g. 
= D t E it 
NEG E py Wier ne 
it 
2 E 1 l 
where IT it at = 
1,t 
with the variable: 
l .= demand for labour 
and as coefficients: 
Kee disembodied technical progress 
As embodied technical progress 
Mo = average productivity of labour in the base year 
8 = number of vintage years of equipment in operation 
A = annual rate of depreciation as a percentage of the 


stock of equipment 


Most important for the working of the model is the fact that 
(i) both the level of production and the rate of growth in 
the preceding years determine the demand for labour, (ii) a 
uniform survival function of equipment for all sectors is- 
assumed and (iii) a linear approximation is applied. 


The total demand for labour each year is calculated as 
follows: 


i 2 
(aux.) l. = Zï Ll. t= 1,...10 


with the variable: 
lt = total demand for labour, year t 
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The total demand for labour cannot exceed labour supply 
minus frictional unemployment. 


(24) 1. < 1s* - uf* 


with the variables: 


ls*t = labour supply in year t 
uf® = level of frictional unemployment 
2.8 Pollution 


Three types of pollution have been included in 
the model. The first group of pollutants consists of 
potentially acidifying substances: sulphur dioxide (SO), 
nitrogen ozides (NO) and ammonia (NH,). The emission of 
these substances is expressed in the unit of acidification, 
mol H*. The second group covers hydrocarbons. Emissions 
are expressed in thousand kilograms. The third group is 
formed by pollution for which control measures have been 
included in the Indicative Multi-Year Environmental Control 
Programme (IMP-M) 1987-1991, an annually updated action plan 
issued by the Dutch Ministry of Housing, Physical Planning 
and the Environment. The content of this category of “other 
pollution" (in total 17 measures) ranges from noise 
abatement to replacement of PCBs and from soil 
decontamination to controlled waste disposal. The emissions 
are expressed in the model in terms of pollution abatement 
costs (in guilders). A distinction is made between emission 
before and after abatement. 

With regard to abatement five sectors are 
distinguished. The first three - sectors 19, 20 and 
21 - concern the abatement of potentially acidifying 
substances in increasing order of cost per unit of 
acidification abated. Sector 22 relates to the abatement of 
hydrocarbons and sector 23 to "other pollution”. 

For the three groups of pollutants, emission 
coefficients have been included in the model enabling the 
pre-abatement level of pollution to be related to production 


We = > 


(25) 


(26) 


3 9 = 


and/or consumption of the goods and services produced in the 
different sectors. In addition emission of pollutants may 
take place independently of the level of production or 
consumption. Thus pre-abatement emission can be defined as 


follows: 
23 18 
S VX ve _ 
ERE EE 7 de Tat Cit KE Dies? 
t= 1,...10 
+ Mt 
with the variables: 
Vk, t = pre-abatement emission of type k in year t 
v*kt = exogenous emission of type k in year t 
and the coefficients: 
Z | 
a = emission coefficient for production 
ae emission coefficient for consumption 
| L 


We will now define the emission level after abatement, i.e. 
after actions are undertaken in sectors 19-23. No additional 
abatement, that is on top of the current abatement can be 
achieved in the first year as the capacity of the specified 
additional abatement sector at the end of the base year is 
zero (see equation (3)). So for year 1 holds 


with the variable: Ê 
Zk,1 = post-abatement emission of pollutant type k 


From year 2 onwards the pollution can differ from the 
pre-abatement emission because of production in the 
abatement sectors. As regards potentially acidifying 
substances three such abatement sectors have been 
distinguished, thereby doing justice to the substantial 
difference in costs between the various control measures. 
The costs per mol H* abated are specified for a number of 


(27) 


` (28) 


= 20 -= 


representative methods in each sector, thus enabling the 
technical, capital, depreciation and labour coefficients as 
well as the import ratio to be calculated. 

Not all annual emissions can be abated. by 
these sectors. The total pre-abatement emission is divided 
into four segments, indicating the technical ceilings of the 
three abatement sectors and an unabateble part. 


Z = 
Vi tři tů Tyt Mog i= 19,...21 
t= 2,...10 
with the coefficients: i 
` 5 
Or = abatement of pollution (for i = 19, ...21 
acidifying substances in mol,H*) per guilder of 
gross output of sector i | 
T; = shares of total emission abateble by sector i 
Post-abatement emission of acidifying substances can now be 
defined as 
21 z 
Ze = Ms - = re t= 2⁄4 0.410 
i=19 


With regard to hydrocarbons and "other pollution" no 
distinction in the cost-effectiveness of different kinds of 
abatement techniques has been made. The possibility of an 
unabateble segment of total pollution however is left open: 


(29) 


(30) ` 


(31) 


= FL se 


Z 
< . PAn 
di teži tó Ti t Vi-20,t i= 22,23 
t= 2.7.4 810 
with the coefficients 
T = share of total emission abateble by sector i 
Z ë . : 
q, = abatement of pollution (for 1.= 22 hydrocarbons in 


grams; for i = 23 “other pollution" in guilders) 

per guilder of gross output of sector i. 
Post-abatement emission of hydrocarbons and “other 
pollution" is defined as gross emission minus abatement: 


= = - t= 2,...10 
Zok =" Vaa Toen Tige “iie i 
ZH _ _ ZH 
Loo,t*z20,t Forse *21,8 
z = V = q ? t= 2,...10 
3,t 3,t 23,t*23,+t 


with the coefficients: 


abatement of pollution per guilder of gross CHEPU 


yé 
SH of sector i 


; = abatement of hydrocarbons (in grams) per guilder of 
/ output of the acid abating sectors i = 19, ...21. 


In equation (30) the activities of the acid-abating sectors 
appear because some of the measures against acidifying 
substances also reduce the emission of hydrocarbons (e.g. 
catalysts in motor vehicles). 


2.9 Goal variables 

An unlimited number of goal variables can be 
formulated: all possible linear combinations of variables in 
the model. In practice we specified nine of them of which 


(32) 


— ZZ =— 


seven were normally used in the optimization process 
(modelling the emission of hydrocarbons ran into data 
problems due to the difficulty of finding an adequate 
summation measure for the different types of hydrocarbons). 
I. Consumption 

Average total consumption, domestically 
produced and imported, is labelled goal variable C. 


io 18 10 18 
i MN 
Cs ta) E Ec +4 E ït c. 


with as variables: 
(e) = consumption by households and government, 
domestically produced 


and, as coefficients: 
MN . . . 
q . import coefficient for non-competing consumer goods 
and services 
qe import coefficient for competing consumer goods and 
services 


2% nsumption rn (F 

Not only can the total volume of consumption 
be an objective, but a balanced development pattern of the 
consumption per sector and between sectors can also be 
desired. A “balanced development pattern” is here made 
operational by defining an additional variable as the 
“lowest” growth rate of consumption in any year in any 
sector. "Lowest" is defined as the deviation from the lower 


limit set, taking into account the distance of this lower 


limit to the upper limit. This is done by expressing the 
deviation as a percentage of the distance. This percentage 
deviation can be calculated for each sector and each year. 
The minimum value of this set of deviations is the variable 
to be maximized. 

For example: assume that the consumption from 
sector 1l in a certain year grows by 2%; the lower limit 
being -2% and the upper +6%. The deviation from the lower 
limit is thus (2-(-2))/(6-(-2)) = 50%. Call this percentage 


(20) * 


(34) 


Fi: The consumption from sector 2 growths by 3%, its lower 
limit being 0% and upper limit 9%. The deviation, expressed 
as a percentage of the distance between lower and upper 

limit, is for sector 2 (3-0)/(9-0) = 33% = F„. Define F as: 


F = min (F, F2) = 33%. 


In maximizing F the consumption from sector 2, being the 
"lowest", is the relevant growth rate. Extending this 
principle to the 18 sectors and ten years prospected, 
equation (20) is reformulated as : 


Eyle, i= 1,...18 i#12 


L 
eren et BCe 7 Yie’ Cio 
1 T 1 


i,t 


with the variable: 
C = consumption by households and government, 
domestically procuced 


and as coefficients: j 
X = lower limit for consumption growth 
y= upper limit for consumption growth 


F is the goal variable to be maximized. 


3. Demand for labour (L ; 

j Average total demand for labour, an indicator 
for employment, is the third goal variable specified. Its 
definition is straightforward. 


with as variable 
l. = demand for labour in year t 
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4. “Highest” unemplo n W 

As with consumption, not only total employment 
but also the course of unemployment in time can be an _ 
objective. Therefore a variable, representing the largest 
deviation from a pre-determined trajectory, is defined as a 


da- >é > 


goal variable to be minimized. 


x = 
(35) Wal - 1. t= 1,...10 
with the variables: 
1, = demand for labour in year t 
1*, = target value for employment in year t 


5. Exports . 
Average total exports E is the fifth goal 


variable. 
L 10 

( 36) E = To 5 e, 
, t=1 


with the variable: 
et= exports in year t 


6. Labour productivity 

The labour productivity of a country is an 
indicator for international competitiveness and wealth. 
Increasing this productivity can form an objective of 
economic policy. To avoid non-linearity this is 
operationalized by formulating a goal variable equal to the 
difference between the index of total value added and that 
of employment at the end of the projection period (average 
over the years 8-10). (Both indices are set at 100 in the 
base year). Value added is defined as gross output minus 
inter-industry sales and imports. | 


(38) 


(39) 


(40) 


= Ge 


( 1 10 ` 23 a 18 . l M 
P = i>i ( ï (1-q, - Ï a, , x. .)}/y 
*t=8 IL a A, eE ist 1,j,t jt [e] 
L 10 . 
Ste 5 Lat 
"tag © T 


with the (exogenous) variables: 


Yo = total value added in the base year 

‘lo = ‘total employment in the base year (thousand. 
man-years) 

7-9. Acidification (Z), Hydrocarbons (H) and "Other 


pollution" (G) 
The average emission of unabated potentially 
acidifying substances, of unabated hydrocarbons and of 


-unabated "other pollution” are the three environmental goal 


variables. 
A 10 
Z = — Ï Z. 
4 10 
H = — Ï z 
10 a 2E 
A 10 
G = — £2 
10 tel 3,t 


with the variable: 


Zkt = nen emission of pollutant k in year 


The model is completed by setting upper and lower Limits on 
all goal variables. 
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NOTES 


(1) Netherlands Scientific Council of Government 
Policy, A Policy-Oriented Survey of the Future; 
Report to the Government no. 25, Staatsuitgeverij, 
The Hague, 1983. 
English summary available. 


Annex 1 SECTORAL BREAKDOWN 


Agriculture (incl. horticulture, forestry and fisheries) 


Mining and gas distribution (oil and natural gas extraction and 
exploration, other mining, natural gas distribution companies) 
Chemicals (basic and endproducts of the chemical industry, 


rubber and plastics processing industry) 

Heavy industry (basic metals, construction materials, 
earthenware, glass and glass products) 

Metal-processing industry (metal products and machinery, 
transport equipment, instruments and optical goods) 
Electrotechnical industry 

Foodstuffs industry (incl. beverages and tobacco) 
Petroleum industry (oil refineries and manufacture of 
petroleum and coal products) 


Other industry (textiles, clothing, footwear and leather goods, 
timber and furniture industry, paper and cardboard, paperware 


and corrugated cardboard, printing and publishing) 
Electricity generation and water supply 

Construction and installation on construction projects 
Trade (wholesale and retail trade plus intermediaries) 
Consumer services (hotels, restaurants, cafés, repair of 
consumer goods, social services, private households with 
Salaried staff) 


. Transport and communications 


Housing 
Business services (incl. 
Health care and education (incl. 


banking and insurance) 
veterinary services) 


. Public administration and defence (civilian and military) 


Additional sectors: 


19. 
20. 
21. 
22. 
23. 


‘Low cost‘ acidification abatement sector 

‘Medium cost‘ acidification abatement sector 

‘High cost' acidification abatement sector 
Hydrocarbons abatement sector (not used in this study) 
Other environmental measures. 


Row nr. 
National 
Accounts 


32, 


34-36 


Standard 
Industrial 
Classification 
(SBI) 1974 


01-03 


12, 19, 402 
ll, 29-31 


32, 33 


22-27 


401, 403 
5 
61-66 


67, 68, 91 


929, 94-99 


71-77 


921-928, 93 


906, 907 


= 728 = 


Annex 2 THE LABOUR DEMAND FUNCTION 


Equation (22) in section 2.7 is derived as follows: 


The demand for labour is a function of production in the 
individual sectors: 


i,t 
(1) Li” 
Mit 
Where : 
Lie = labour demand in year t in sector i, 
Xie = value of production in sector i in year t. 


H = labour prođuctivity of the capital goods 
i,x : Ê ges 
used in the production of sector i in year t. 


As the capital goods used have different productivities 
according to their age and as we want to introduce technical 
progress, we can rewrite (1) as follows: 


À (= t, pe rl 
z > 
(2) Bens al 
Fles T ri 
Where 
Xp; = the value of production in year t in sector i, 
d 1 


produced with capital goods ay years old. 
Led = labour productivity of capital goods used in 
production in sector i, in year t, with capital 
goods T years old. 
6 = number of vintage years of capital goods used 
in the production process. 


As we want to take into account embodied and disembodied 
technical progress influencing labour productivity, we 
specify the following equations that show how labour 
productivity increases or decreases over time: 


(3) Me N oe $" 
(4) Me = Mo < Gt Dee 3 


GIN Moy = Mos Aden AQ I oe 

Rewriting (5) with y~ for the years in the past (and t 
for the years in the future) and combining (3),(4) and 
(5), we obtain the following equations, where for reasons 
of clarity we omit index i: 


(6) tea Mo .( l+ Àp yt ( 1+ p E tT 
Where j | j 
Mo = labour productivity in the base year preceding 
the prospected period. 
AE and À p = coefficients representing embodied and 
disembodied technical progress. 


Substituting equation (6) in equation (2) results in: 


9 X 
tT 
(7) L t = 2. 
T =1 Mo . (1+ A yt Er Àe yt-T 


Multiplying the right-hand side by 1 / (1+ Ap), 

rewriting it for period t-l and then subtracting this new 
expression from both the left- and right-hand side of 
equation (7), we obtain: 


L X 
(8) L t= = + t, = 
l+ t t-1 
D Ela (1+ À p) (1+ À g) 
Steine 
t t-1- 8 
j 
Le . (l+ Ap ) (1+ Ar ) 
because X t,2 = X Ent 
and XL 3 = Ke-1,2 , etc. 
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and the expressions between the first and the last ones 
cancel each other out. 


Equation (8) shows that the labour demand in year t is 
reduced by disembodied technical progress (first term on 
the right-hand side), but may expand because of (new) 
production made possible by access to new capital goods 

( T =l, second term on the right-hand side), and again is 
reduced by dismantlement of capital goods (from the oldest 
year, T = 8 , the third term on the right-hand side) 
with lower labour productivity per output. 


Rewriting equation (8) in order to relate the demand for 
labour directly to the total level of production, and not 
to various levels of production made with different 
vintages of capital goods, we must introduce x, 
representing the percentage of capital goods set aside 
annually. 

Therefore we define: 


(9) Xt,1 = Â Xt + x, Xt-1 
Equation (8) can be rewritten as: 
D Nort Ap) Cs ^ p) 
X t-1, 6 


Hog (lt Ap E Ge Api? 


We note that : 


(11) Kil, 0 = % , x 


or, in words: the oldest part of the capital goods still 


9 


in use (from year ) equals that part that is set aside 


(thus & ). Using (11), equation (10) becomes: 


Le 4 AX + La ., X. 1 


N g(t À 0) E14 Age 


< `+ *#t-1 


Ho (+ 2 DEC À peo? 


Equation (12) shows the labour demand in period t as a 
function of the labour demand in period (t-1), and also as 
a function of the production increase and of production in 
the previous period. As production increases, so will 
labour demand. It can be rewritten as follows: 


bea 
(13) L = + 
l+ À D 
Ə 
AX + x Kef 1-(1+ A E) ) 


Mo ot A pede Ag I 


A practical problem is that the labour productivity of 
capital goods produced in the base year (the term 1 9) 

is not known and we want to write it in terms of the 
average labour productivity in the base year. This average 
equals: 


(14) Mo = = 
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Substituting the Mo as presented in (14) in equation 
(13), we obtain: 


L 
(15) Ly = t-1 r: 
9 i>] 2 
AX.( (1+Àp) -1) - SR CX 1+ Àg) -1 ) 
+ 


no Be Ag (1+ À D j". (1+ À E -l+ ° 


This expression equals: 


L 
(16) bys ECE =< ANER + 
1+ À p 
(=) 
. DE Ft Re oe Cl Se Ci Xg) ) 
= t t l- g 
Ho 9 (l+ AD). (lt À œ) `+ - $ 
-g 
l 
Where g = ee 
l + p` E 


Equation (16) is the same as equation (22) presented in 
the main description of the model. 
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Annex 3 A COMPARISON BETWEEN LINEAR PROGRAMMING IN 
GENERAL AND REVISED SIMPLEX 


In general 


Linear programming consists of solving the following 
problem: 

maximize a variable Z which is a linear function of n 
variables x(1), x(2),...,x(n).. 


(1) Z = c(1).x(1) + c(2).x(2).+ ... +c(n).x(n) 
(Z is the objective) 
while at the same time fulfilling a set of conditions: 


(2) a(1,1).x(1) + a(1,2).x(2) + ...+ a(l,n).x(n) < b(1) 
a(2,1).x(1) + a(2,2).x(2) + ...+ a(2,n).x(n) < b(2) 


a(m,1).x(1) + a(m,2).x(2) + ...+a(m,n).x(n) b (m) 


Ja 


These are the m restrictions. 


In matrix notation: 


Where Z is a variable (scalar), 
C is a row vector of coefficients c with length n; 
X is a vector of variables with length n; 
A is an mxn matrix of coefficients; 
B is a vector of fixed values with length m. 


First we introduce the m slack variables x(n+l),...,x(n+m) 
in order to transform the inequalities (2) into equations: 


- (3) A.X+#+tI. Xe = B 


= 34 = 


Where I is an m x m identity matrix, 
Xs is a vector of length m. 


The system of equations (3) can be rewritten as follows: 


(4) (A Í I) En á ` `> = 
Xs | 


This is the set of equations that has to be solved. With a 
set of X and Xs that solves this system, one can 
calculate the value of Z, the objective. One can then see 
if this Z is optimal. If not, one has to find another set 
of X and Ke that solves (4) and calculate Z again. This 
procedure is continued until the value of Z cannot be 
improved by a different set of X and Ke 

This search process is equivalent to the following 
iterative algorithm: 


Find one set of X and Ks 
that solves (4) 


Calculate Z, the 


objective 


Is Z optimal? yes STOP 
no 


Find another set of 
X and Xs that solves 
the system of equations (4) 


In geometrical terms this search method means that the 
linear restrictions define and enclose a space for the 
variable Z = sum of c(i). x(i). In the process of 


searching for a solution, the contours of this space are 
explored. 
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The larger the number of X variables, i.e. the larger n 
is, the greater the dimensions of the space that has to be 
explored. Therefore, the larger n+ m is, the number of 
variables, the larger the number of iterations that can be 
expected to be necessary. On the other hand, the larger m 
is, the number of restrictions that have to be met at each 
iteration, the longer the calculations take during each 
iteration to find a solution that satisfies. 


In algebraic terms, the solving of a system of equations 
is equivalent to inverting a (in this case) mxm matrix and 
requires approximately m? additions and multiplications. 
This is how the computations are done in the Simplex 
method. In our model we use the Revised Simplex method. In 
the Revised Simplex method the calculations are shortened 
and simplified as much as possible as compared to the 
Simplex method. In the Revised method this is done by 
defining and distinguishing between non-basic variables 
(the ones that are equal to zero) and the basic variables 
(the ones that are not equal to zero). This simplifies the 
calculations because we will not include all n+m ( X and 
XS) variables in the calculations but only the m basic 
ones. We define: j 


1 


‚BAS = ( A | I) but only with the m 
columns corresponding to the basic variables. 


Also: 


XB = ( X i Xs )` but only with 
the m basic variables, 


+ 
t 


CB = ( C ! 0 ) but only with the 
m corresponding basic c's. 


The problem can be reformulated as follows: maximize (5) 
under the restrictions (6) 


(5) Z = CB. XB 

(6) BAS . XB = B 

Equation (5) replaces (1) and (6) replaces (3). 
Equation (6) can be rewritten as 


(7) XB = BAS Í . B 


This is the lution for the X basic variables after each 
iteration. j 
Equations (5) and (7) combined 


(8) Z = CB . BAS . B 


This is th lution for Z after each iteration. 


Again the above equations can be rewritten. Taking 
equations (7) and (8) together: 


Z \ B CB BAS lÈ . B 
(23 XB BAS l. B 
OR 
3 ! =| , 
z \ _ 1 ' CB. BAS ) 0 
: 1___1_ CB. BAS) 
(10)\ xp 0 i pas Í B 


Equations (1) and (3) can be put combined as: 


0 1 -C - 0 Z 
(11) = . X 
B 0 A I X 


s 


=37 = 


Replacing the ( 0 :B ) vector in (10) with its value in 
(11): 


Z L Í CB.BAS”Ì a -c | cB. Bas !| /z 

(12) Ss de < = > = >> - > > < < => = ee x 
i -1 ! = | 

XB o | BAS LA ' BAS xs 


Using equation (9) for the left-hand side of equation (12): 


cB. BAS? . B 
(13) F. oe anaemia E idem 
as 
gi RHS (12) 
BAS l. B 


This is the form the set of equations has after every 
iteration. Note that this equation uses the same A and B 
each iteration, and CB and Bas t only change a little 
after each iteration. This is therefore a handy algebraic 
form to use in calculations. It can be shown that each new 
BAST} 
computation, as after each iteration only one columm 


can be deduced from the former, by simple 


vector of this m size matrix changes. This "Revised" 
Simplex method is therefore more efficient than the 
Simplex method where each iteration involves calculations 
very similar in length and size to calculating the inverse 


of a new mxm size matrix. 


Annex 4 LIST 
(Éor 


All variables 
prices unless 


uf* = 


Vk,t = 


Wi,t = 
Xi,t = 
Yo = 


Zk,t = 
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OF VARIABLES INCLUDED IN THE MODEL 
mođe of notation see 2.1) 


are expressed in million guilders, constant 
otherwise stated. 


consumption by households and government of 
goods and services produced by the domestic 
industry i in year t. 


exports of goods and services produced by 
industry i in year t 


total exports of goods and services in year t 


gross investment in fixed assets produced by 
industry iin year t 


exogenous investment produced by industry i 


demand for labour in industry i in year t 
(thousand man-years, constant working hours) 


total demand for labour in year t (thousand 
man-years, constant working hours) 


labour supply, in year t (thousand man-years, 
constant working hours) 


target value for employment in year t 
(thousand man-years, constant working hours) 


imports of goods and services of sector j in 
year t (excl. non-competing imports of 
consumer goods and services) 

total imports of goods and services in year t 


interest margin of banks (not equal to zero 
only for sector 16: Business services) 


frictional unemployment (thousand man-years, 
constant working hours) 


pre-abatement emission of pollutant k in year 
t (106 mol Ht for k=l : 103 kg for k=2; 

10° guilders for k=3) 

production capacity of industry i in year t 
gross output produced by industry i in year t 


total gross value added in base year 


post-abatement emission of pollutant k in year 
t (units: see vk,t above) 


Goal variables 
C = 
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average total consumption, domestically 
produced and imported (Cl and CU lower and 
upper limits) 


indicator for a balanced consumption pattern 
(unitary dimension; FL and FU lower are -- = 
upper limits) 


average demand for labour (thousand man-years; 
LL and LU lower and upper limits) 


largest deviation in any year from a 
predetermined target trajectory of employment 
(thousand man-years; WL and WU lower and 
upper limits) 


average exports (EL and EU lower and upper 
limits) 


indicator for rise in labour productivity 
during the projection period (unitary 
dimension; PL and PU lower and upper 
limits) l 


average unabated emission of potentially 
acidifying substances (106mol H+; zl and 
zU lower and upper limits) . 


average unabated emission of hydrocarbons 
(10kg; HE and HU lower and upper limits) 


average unabated emission of “other pollution" 
(GL and GU lower and upper limits) 


Annex 5 


a. … 
1,j,t 


d. . 
i,j,t 


1,¢€ 7t 
D E 
À 
i,t’ ^ , t 
u. 
1,0 
8 
a 
T, 
1,t 
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LIST OF COEFFICIENTS INCLUDED IN THE MODEL 
(for mode of notation see 2.1) 


technical coefficient: sales of industry i 
needed to produce one unit of output by 
industry j in year t 


depreciation coefficient: sales of investment 
goods produced by industry i needed to 
maintain one unit of production capacity in 
industry j in year t 


capital coefficient: sales of investment goods 
produced by industry i needed to expand the 
production capacity in industry j by one unit 
in year t. 


minimum average degree of utilization of 
capital in industry i 


maximum permitted annual expansion of 
production capacity in sector i in year t (as 
% of capacity in the previous year for 
i=l,...18; in millions of guilders. for 
1219,... 23) 


upper or lower limit of change in exports of 
sector i in year t, expressed as a percentage 
of exports in the base year 


upper or lower limit of the balance of trade 
in year t (millions of guilders) 


upper or lower limit of the average balance of 
trade during the projection period (millions 
of guilders) 


upper or lower limit of change in consumption 

of goods and services produced in the domestic 
industry i in year t. Expressed as a ratio to 

domestically produced consumption of industry 

i in the base year 


disembodied or embodied technical progress in 
industry i in year t 


average Labour productivity in inlet i in 
the base year 


number of vintage years of equipment in 
operation 


annual rate of depreciation of equipment as 
percentage of the stock of equipment 


share of total gross emission of pollutants 
abateble by sector i 


ZH 
{i,t 
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import coefficient of expansion investment, 
expressed as a ratio to total domestically 
produced expansion investment in the recipient 
industry j in year t 


import coefficient of replacement investment, 
expressed as a ratio to total domestically 
produced replacement investment in the 
recipient industry j in year t | 


trade margin on the exports of goods 


import coefficient of raw materials and 
semi-manufactures as a ratio to domestic gross 
output in industry j in year t 


import coefficient consumption goods and 
services (excl. non-competing consumer goods 
and services), expressed as a ratio to the 
domestically produced consumption of the 
competing industry j in year t 


import coefficient non-competing consumer 
goods and services, expressed as a ratio to 
the total yearly domestically produced 
consumption 


trade margin on consumer goods 


emission coefficient production in industry i 
in year t of pollutant k (mol Ht per guilder 
for k=l; grams per guilder voor k = 2; unitary 
dimension for k=3) 


emission coefficient consumption of. products 
produced by industry i in year t of pollutant 
k (units: see qYX above) 


abatement of pollution per guilder gross 
output of sector i in year t (mol Ht per 
guilder for i = 19,...21; grams per guilder 
for i=22; unitary dimension for i=23) 


abatement of hydrocarbons per guilder output 
of acid-abating sector i (grams per guilder) 
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Annex 6 SEQUENCE OF EQUATIONS IN THE COMPUTER LISTING (23-6-87) 
Label Number Corresponding Type Description 
of rank no. in of 


equat. documentation equat. 


AAAL#.. 7 - E Definition of goal variables 
AAAO#.. 7 (32b), (40c), LE, - Lower and upper limits of 
(34b), (35c), GE goal variables 
(36b), (37b), 
. (38c) 
AGRE#.. ` 7 (32c), (40b), LE, Lower and upper limits of 
. (34c), (35b), GE goal variables 
(36c), (37c), 
(38b) | 
X...#.. 180 (1) - E Leontief equations 
W...#.. 230 (2) LE Production cannot exceed 
capacity 
J...#.. 230 (4), (5), (6) LE, E Maximum expansion of 
capacity 
I...#.. 200 (7) LE Maximum contraction 
of capacity 
I...#99 23 (3) GE Minimum utilization 
of capital 
GAAA#99 1 (8) GE Investment in last year 
E...#.. 170 (9) LE Upper limit for. exports 
EALL#.. 10 (12) E Total annual exports 
EALL#99 1 - E Total exports 
F...#.. 160 (10) GE Lower limit for exports 
FALL#.. 10 - E Total annual consumption 
FALL#99 1 - E Total consumption 
EHAN#.. 10 (11) E Export trade margins 
M...#.. 230 (13) E Imports 
MALL#. . 10 (14) E Total annual imports 
MALL#99 1 - E Total imports 
PAAA#99 1 (17) LE Upper limit for average 
balance of trade 
PBBB#99 1 (18) GE Lower limit for average 
balance of trade 
HAAA#.. 10 (15) LE Upper limit for annual 
balance of trade 
HBBB#.. 10 (16) GE Lower limit for annual 
balance of trade 
C...#.. 170 (19) LE Upper limit for 
consumption 
D...#.. 170 (20) GE Lower limit for 
consumption 
A000#99 1 - E Scaling of goal variable F 
FHAN#.. 10 (21) E Consumption trade margins 
L...#.. 230 (22) E Demand for labour 
LALL#.. 10 (23) E Total annual demand for 
labour 
LALL#99 1 - E Total demand for labour 
OAAA#.. 10 (24) LE Employment cannot exceed 
labour supply 
V...#.. 30 (25) E Pre-abatement pollution 


oe ee 


ZAAA#.. 
ZBBB#H.. 
ZCCC#.. 
COOO#99 
LOOO#99 
WOOO#.. 
EOOO#9 9 
ZDDD#.. 


ZEEE#.. 


HALL#99 
HOOO#9 9 


HAAA#99 
GALL#99 
GOOO#9 9 
GBBB#99 


TALL#.. 
TALL#99 
POOO#99 


ZALL#99 
ZOOO#99 


u) 
oo 


OHHH 000 
oo 


wo 
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(26), (28), 
(30), (31) 
(27) 

(27) 

(27) 

(32a) 

(34a) 

(35a) 

(36a) 

(29) 


(29) 


(39a) 


(39c) 


(40a) 


t tj 


tm AM 


Post-abatement pollution 


Upper limit acid-abatement 


Goal variable consumption 
Goal variable employment 
Goal variable unemployment 
Goal variable exports 
Upper limit for hydrocarbon 
abatement j 
Upper limit for 
"other-poll." abatement 
Unabated hydrocarbons 

Goal variable 

hydrocarbons 

Upper limit goal var. 
hydrocarbons 

Unabated "other pollution" 
Goal variable "other 
pollution" 

Upper limit goal "other 
poll." 

Total annual value added 
Total value added 

Goal variable 

productivity 

Unabated acid 

Goal variable acid 


